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The hadronic interaction model BBL implements the ideas of gluon saturation due to large densities. When 
approaching the black body limit at high energies, leading partons acquire large transverse momenta which 
breaks up their coherence. This leads to a suppression of forward scattering, and is therefore important for air 
showers. We discuss some general aspects of this new approach and their influence on air shower properties 
as seen by fluorescence and surface detectors: The position of the shower maximum is reduced due to stronger 
absorption in the atmosphere. The lateral distribution functions become flatter for the same reason. Muons are 
produced abundantly due to high multiplicities in the mid-rapidity region. The response of water Cherenkov 
detectors and comparisons to other interaction models are shown. 


1. Introduction 

In this paper, we discuss the high energy limit of hadron nucleus scattering and their influence on air shower 
properties, as described in Ref. [H. When approaching highest energies, we expect the differential scattering 
amplitude to become close to unity. The partons acquire a large transverse momentum of the order of the 
saturation scale, which leads to an steeper spectrum of forward scattered particles, the most important phase 
space region for air shower properties. 

All air showers in this paper have been computed with the Seneca model |2| using UrQMD 1.3.1 O as low- 
energy model below 100 GeV. 
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Figure 1. Left panel: Xmax as a function of energy. Right panel: Xmax of models using inelastic cross-sections from 
QGSjetOl, i.e. differences are due to forward scattering only. 
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Figure 2. Left panel: Inelasticity of the three models as a function of primary energy. Right panel: Inelasticity Kq.qi, 
evaluated in the forward phase space region 0.01 < a;F < 1. 


2. High energy hadron nucleus scattering 

When approaching highest energies, higher twist corrections in hadron nucleus collisions become increasingly 
important. Attempts to account for this are the implementation of an energy-dependent pt cutoff m for hard 
scattering or the resummation of enhanced pomeron diagrams in an efficient manner (5). 

Our approach is to consider the black disk limit (or black body limit - BBL) at high gluon densities within the 
Color Glass Condensate (CGC) approach (6). where the interaction probability is close to unity. Within this 
model, the scattering amplitude of a quark on a dense gluon field is resummed to all orders. This can be done 
since the typical scale of the gluon density, the so-called saturation momentum Qg is larger than Aqcd and 
weak coupling methods are applicable. 

The Monte Carlo interaction model, which incorporates the physics discussed here has been introduced in 
Ref.lH and discussed in more detail in Ref. |2|. The most important feature is the suppression of forward 
scattering. As a consequence of energy conservation we also find an increase of particle production at lower 
rapidities. Although less important for the longitudinal profile of an air shower, this results in an increased 
muon production and is therefore important to be considered when deducing composition from this observable. 

Since this approach is only valid when the saturation momentum is large, we use this model only in the high 
energy/high density limit. For low energies and peripheral interactions we use the standard pQCD Monte Carlo 
model Sibyll |4). 


3. Longitudinal Profile 

Results for the mean Wmax are shown in Fig. ^ The longitudinal profile is dominated by mainly two observ¬ 
ables, the forward scattering, often described by the inelasticity (K = 1 — of the fastest particle) and the 
inelastic cross section which determines the mean free path. Since we have used the inelastic cross sections 
from Sibyll, differences to this model are due to the reduced forward scattering solely, whereas the difference 
to QGSjetOl II comes from a combination of both different cross sections and inelasticity. 

Also shown are the results for iron induced showers computed with the BBL model by simple superposition 
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Figure 3. Left: Rapidity of charged particles for p+Air at 10^^ eV. Right: Electron/positron and muon numbers in the 
intermediate energy range. 


of nucleon-air collisions, just as in Sibyll. This should give a good estimate, since the the running coupling 
evolution is less sensitive to initial conditions, i.e. the density of the projectile. The full implementation of 
nucleus-nucleus collisions with consideration of projectile saturation scale will be left for the future. 

Fig. El shows the inelasticity of the models, defined as mean 1 — xf of the fastest particle. The right panel 
shows iTo.oi = 1 — Jqqi xpAn/dxpdxF, which corresponds to the inelasticity averaged over a larger region 
in forward phase space, not only the fastest particle. 

The right panel in Fig. [Q shows Xmax for Sibyll and BBL when one applies inelastic cross-sections from 
QGSjetOl which are somewhat smaller, see e.g. O. At largest energies the Xmax of Sibyll and BBL get 
shifted by ca. 15 g/cm^ and reach 50 and 30 g/cm^ larger depth than QGSjet, correspondingly. Assuming the 
same inelastic cross-sections, BBL has a larger X^ax than QGSjetOl, despite a larger inelasticity of the former. 
But as can be seen in Fig. |2] QGSjetOl has a larger Xo.oi, as a consequence of a rather flat x-distribution at 
these energies. Clearly, the inelasticity alone is not sufficient to characterize forward scattering. In Ref. Q is 
was shown that X„iax is sensitve to the forward region down to xf ~ 10“^ — 10“^. Therefore the Xo.oi, or 
even Xo.ooi variable is more suitable to describe this specific range. 


4. Muon- and electron numbers 

The major consequence of the implementation of the black body limit is the suppression of forward scattering. 
However, this also leads to an enhancement of the multiplicity in the mid-rapidity region, see Fig. 0 This 
has almost no influence on the longitudinal profile, but increases the muon numbers as decay products of low- 
energy particles. A plot of mean electron/positron and muon numbers is shown in the right panel of Fig. 0 
The BBL model seems to interpolate somewhat between the results of Sibyll and QGSjetOl, with the transition 
happening at about 10^^ eV. Since the primary energy is shared between the nucleons of the projectile nucleus, 
this interpolation also happens as a function of mass number A at correspondingly higher energy. This can be 
seen for the carbon projectile in the knee region. Whether this model can help to resolve known problems of 
models in KASCADE, is subject to a detailed comparison. 
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Figure 4. Lateral distribution functions of different particles and of the total signal in VEM units. Results are shown as a 
ratio to the LDFs from Sibyll. The low energy model (E < 100 GeV) is UrQMD. 


5. Lateral Distribution Function 

In Fig. 12 we show the lateral distribution functions for 10^® eV vertical proton showers at Auger altitude. 
Ratios to the LDFs of Sibyll are plotted. Besides from a small decrease of the slope for electrons and photons, 
most noticeable is the increased muon number. QGSjetOl shows the same qualitative behavior as BBL. Since 
muons dominate the VEM signal at large distances, the combined LDF in VEM units shows a signihcantly 
flatter slope. 
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